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Abstract 

We measured the polarized and azimuthal dependencies of optical second har- 
monics generation (SHG) on polished surfaces of ZnSe(lOO) single crystal sur- 
face in air, using a fundamental wavelength of 1.06/um. By considering both, 
the bulk- and surface- optical nonlinearities within the electric dipole approx- 
imation, we analysed the data for all four combination of p- and s-polarized 
incidence and output. The measurement using Sin — Sout is thereby partic- 
ularly useful in the determining of the symmetry ZnSe(100)/contamination 
layer-interface, which would lower the effective symmetry of ZnSe(lOO) from 
Civ to C2v The analysis of p-incident and p-output configuration allows us 
to distinguish the [Oil]- and [Oil]- directions. 

I. INTRODUCTION 

Optical Second Harmonic Generation(SHG) has matured to be an effective surface- and 
interface-sensitive technique for probing a large variety of surface- and interface-properties. 
This technique is a non-destructive. It can be used for in-situ measurements in almost 
any environment, provided that we have optical access to the sample [jl],?,?]. The method 
bases on the idea that generally a surface and a bulk have different structural symmetry 
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0,?,?,?]. Materials with inversion symmetry, such as Si- and Ge- single crystals have no 
bulk electric-dipole contribution. Electric quadrupole and magnetic dipole contributions 
originate the leading source in SHG from the bulk. However, a surface or an interface 
breaks the inversion symmetry and produces an electric-dipole contribution. Generally the 
SHG-intensity from a surface or an interface can be compared to be much greater than 
that from the bulk, depending on the material, the photon energy, and the geometry of 
the experiment. On the other hand, noncentrosymmetric materials like crystals with the 
zincblende structure as, e.g. the compound semiconductors GaAs and ZnSe, have intense 
bulk electric-dipole contributions to the SHG, thus hampering surface observations. Stehlin 
et al 1^ illustrated possible combinations of the polarization of incident and SHG output 
light and of crystal azimuths for which only the surface contribution may be deduced for 
typical low index surface, e.g. (001), (110), and (111) and used SHG as effective surface 
probe with submonolayer sensitivity to monitor the adsorption of Sn on GaAs(OOl) surface. 
Later, S. R. Armstrong et al [0 investigated by means of SHG a GaAs(100)/air interface. 
Similar experiments performed by C. Yamada and T. Kimura |^,?] firstly observed the 
twofold rotational angledependence (rotation anisotropy) of SHG in the reflected light from 
a well prepared and by RHEED characterised noncentrosymmetric single crystal surface 
of GaAs(lOO). They interpreted it as the interference of surface SHG and dipole-allowed 
bulk SHG having fourfold symmetry results in twofold anisotropy, whereby the degree of 
rotational anisotropy monitors the surface reconstruction. Galeckas et al showed that 
SHG in reflection is a sensitive and practical method to invesigate the departure from perfect 
crystalline ofder at the surface by studying the SHG rotational anisotropy of crystalline, 
polycrystalline and amorphous silicon carbide surfaces. Bottomley et al determined to 



within ±0.1° the orientation of vicinal GaAs(OOl) and Si(lll) single crystal wafers using 
second and third harmonic generation. M. Takebayash et al [|12| demonstrated the SHG- 
measurement of the tilt angle 9 of the crystallographic axis of a vicinal GaAs(lOO) wafer 
towards a direction ^ that characterises step direction and step height and forms an angle 
ip with respect to [100]-axis and found that the s-incident polarization is useful in the 



determination of this tilt angle. C. Jordan et al used the azimuthal rotational anisotropy 



of the SHG for in-situ finger print characterisation of various polytypes of SiC and compared 
with X-ray diffraction (XRD) and identify defect regions in big samples by observating the 
spatially resolved dependence of the SH-intensity. 

In recent years, there has been big interest in wide-band gap II- VI compound semi- 
conductors. In particular, ZnSe is being actively studied due to their electroluminescent 
properties in blue-light-emitting laser diodes [?]. These devices are fabricated by epitax- 
ial growth on GaAs substrates. The replacement of GaAs substrates by ZnSe ones may 
decrease the concerntration of defects and increase the lifetime of the device. This can 
be one possible way to improve ZnSe-based laser characteristics. An recent observation in 
ZnSe/ZnCdSe laser diodes grown onto ZnSe substrates indicates that the growth rate of the 



defects is much smaller than for layers grown onto GaAs substrates |T^. For the successful 
growth of II- VI epilayers onto ZnSe layers is very importment to have well ordered substrate 
surfaces. However, the superiority of ZnSe substrates over GaAs one has not been realised 
yet, mainly because it is difficult to obtain high-quality ZnSe substrates. 

SHG has been successfully apphed to the study of semiconductor surfaces semiconduc- 
tor/semiconductor heterostructures and interfaces, and semiconductor-oxide interfaces |I| 

and provides a powerful tool for understanding the energetics of various defect formations 
and surface reconstructions. Although these experiments have been performed mostly for 
GaAs surface or SiC surface similar symmetric studies were not extend to 



many other non-centrosymmetric materials, such as ZnSe. 

In this paper, we have measured the SH intensity from ZnSe (100) as a function of the light 
polarization and as a function of the azimuthal angle. We have also calculated the intensity 
of SHG from the surface electric dipoles and the the second-order nonlinear polarization 
arising from the bulk response. We use optical SHG to determine the macroscopic symmetry 
properties of the ZnSe(100)/air interface. The remainder of the paper is organized as follows. 
In Sec. II, the theory for SHG in reflection from the bulk and surface of nonsymmetric 
crystals is briefly introduced. Then the experimental apparatus is described in Sec. HI, and 



finally the results are analyzed in terms of symmetry. 



II. THEORY 



A. General 



We can express the induced nonlinear surface polarization at an interface in the following 
manner: 

Pe//(2CU) = P,^rf + PbulkF{uj)F{2uj)L,ff, (1) 

where F{uj) and F{2uj) are the Fresnel factors for the incident input and output fields, and 
-^e// = (W + 2w)~^ is the effective phase- matching distance of the substrate with W = fc^j^z 
and 2w = k2u).z representing the z component of the wave vectors of the fundamental and 
SH light, resepectively, Pswr/ is the surface nonlinear polarization of the interfacial layer, 
and Pbuik is the bulk nonlinear polarization in medium. 

The bulk nonlinear polarization can generally be expressed by a series of multipole terms: 

m ■■ 

Pbuik{2uj) = P^^\2uj) - V ■ Q^^\2uj) + X M(2)(2cu), (2) 

where P^^-*, Q*^^\and M*^^) describe the electric-dipole polarization, electric-quadrupole po- 
larization and magnetic-dipole polarization respectively. Up to the first derivative in Pe//; 
the three polarization sources are the following: 

P(2)(2c.) = : Eiu)Eiu;) + : Eiu)VE{u;), (3) 
Q^'\2co) = x'^:Eico)Eico), (4) 

M'^^\2uj) = x^' ■.E{uj)E{uj). (5) 

We neglect M and higher-order multipole in our discussion here, because these contributions 
are thought to be orders of magnitude smaller than the electric dipolar contribution. The 



second terms in eq. (0) and (|^) arises from the gradient of the field E{uj). Since the 
penetration depth in crystals (~100A) is much smaller than the spatial variation of the 
field (~10,000A for optical frequencies), this bulk contribution is small under nonresonant 



conditions. [p7 



The surface contribution of the SH field arises from two effects. First, at the interface 
between the dissimilar media, inversion symmetry is broken for centrosymmetric crystals 
and symmetry is changed for noncentrosymmetric crystals, and so a dipolar contribution to 
SHG can exist. In addition, there is a discontinuity in the fundamental electric field normal 
to the surface. This can generate a sizable SH contribution through higher-order multipole 
terms. So the surface contribution to the nonlinear polarization can be expressed as: 

Pi2/(2^) = X"" ■■ E{u)E{u) + : E{u)VE{u) - : VE{u)E{u). (6) 

Since the tangential components of the incident electric field are continuous across the 
interface, all terms in eq. (^) involving a gradient that can be neglected when describing an 
in-plane surface response and only the first term should be considered: 

^2/(2^) = x'':^M^M- (7) 

If the nonlinear susceptibilities are all defined in terms of the input field E(ci;), in the 
interfacial layer by ^{oj) = x*'^^ • E{uj)E{uj), then from (|^), the corresponding effective 
surface nonlinear susceptibility is given by 

X?A = xfif + X?],,L,ffF{u)F{2uj), (8) 

where xlurf ^^'^ xluik surface and bulk nonlinear susceptibility tensors, respectively. 

(2) 

The tensorial properties of xl/f can be exploited as: if the SH intensity is recorded as a func- 
tion of azimuthal angle of rotation, the variation in intensity reflects the overall symmetry 
of the surface of interface. 



A tensorial expression for the second-order polarization can be written as: Wl 
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Xxxx X,xyy X,xzz X,xyz X,xxz X,xxy 

v v V v v V 
Ayxx A.yyy A.yzz A.yyz A.yxz A.yxy 



(9) 



E^{u)Ey{u) 

Ey{uj)Ey{uj) 

E,{u)E,{u) 
2Ey{uj)E,{uj) 
2E,{uj)E,{uj) 
y 2E,{uj)Ey{uj) j 

This is the general expression describing the interaction of the two EM driving fields being 
coupled throuth a dyadic product. 



B. Bulk contribution 

We now define a new set of coordinates (x', y', z') for each of the three crystal orientations 
such that z axis is perpendicular to each crystal face. For the (111) crystal face, we therefore 
have, in terms of the standard crystal axes (x, y , z) 

'^x'^ U\ /2/^/6 -1/a/6 -l/x/6\ /^x^ 

^1 y = 1/^/2 y , (10) 

^z j 1/^3 1/^3 / \Zy 

where the new x' axis is projected on to the original crystal x axis in the plane of the 
crystal surface, and (x, y,z) represents the principal-axis system of the crystal ([100], [010], 
[001]). For the (100) face we simply choose the x axis to he normal to the surface, the 
transformation matrix can be expressed by 

^0 1/^2 -1/v^^ 
Ri^ l/^/2 1/V2 , (11) 
1 



and for the (110) face crystal 



^ -1/^/2 l/s/2 0^ 

1 

l/x/2 l/x/2 ) 
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We can define a set of unit vectors for tlie incident liglit beam s", k,and z sucli tliat "s and 
k lie on the crystal face, perpendicular and parallel to the plane of incidence, respectively, 
and z = z' is the surface normal, 

/ . 









k 


= Ro 


y 









sm ( 

COS( 





y 



(12) 



- cos 
sin 

i; 

Here (p is the angle between k and x, the azimuthal angle. 

In many of the experiments, the single surface is rotated about its azimuthal angle (f) and 
the SH response is analyzed with respect to the beam coordinates. Therefore, one needs to 
transform x^"^^ from crystal coordinates into beam coordinates, as shown in figure 1, with 
the appropriate transformation operations. The transformation rule for third rank tensors 
is the following: 



(13) 



l,m,n 



where R = RqRi is the transformation operator from the crystal (x, y,z) to beam ('s, k,z) 
coordinates. 

Wide-gap zinc-blende II- VI semiconductors with symmetry Td(43m) are optically 
isotropic, but do not possess a center of inversion. The bulk second-order susceptibility 
tensor has only one component, Xxyz = Xyzx = Xzxy = d. Using the rule in Eq. (p!3D , this 
tensor for the (100) face in beam coordinates can be written as 



~sm{2(j))d cos(2(/))d 

-cos(20)d -sin(20)d 

cos(20)rf -cos(20)d -sin(20)ci J 



(14) 



and for (111) face we have 
/ 

XbuM = 



— -y/ 3 sm{3(f))d 
I cos(30)(i 
d/V^ 



sm{3(j))d -d/V3 

„ cos (30)d -d/VS 



-d/^/3 2rf/V3 
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2 cos(30)rf ^ 



3 sin(30)(i 




(15) 



The components of the fundamental field in the medium expressed in the beam coordinate 
axes are 

Ek = fcEptp, (16) 
Ez fsEptp, 

where /c,s are the Fresnel factors and tg^p are the linear transmission coefficients for the 
fundamental field, given by 

r _ sing r f 2,1/2 . _ 2cosg , _ 2cosg 

n{2uy ^ ' ' cos^ + n/,' " ncosd + f,' 

with fc taken such that Im fc > 0, and Re /c > if Im fc = 0; n{uj) is the complex 
refractive index of the medium in which the beam is propagating. 6 is the angle of the beam 
propagation in this medium. For the SH field, the coressponding Fg^ Fc, Tg, Tp are similar 
with refractive index A^(2ci;) of SH light. 

Through combining Eq. (P) and Eqs. ([T^)-(|l7|), one obtains the second-order polarization 
from the (100) face crystal, 

pbuik ^ _2rf[sin(20)/J,t2^2 _ cos{2<p) UptgEpEg], 

pbuik ^ -2d[cos{2<p)fJgtlEl + sm{2(j))fgtptgEpEg], (18) 
Pf''^ = d[cos{2<p)tlE^^ - cos{2<f))fXE'p - 2sm{2<f))UptgEpEg], 

and for the (111) crystal, 

= -^d[sm{3<P)tlE'^ - sm{3<jy)fXE'p + 2(cos(30)/e + V2fg)tptgEpEg], 
Pfc"'' = -^d[t' cos{3(P)E^g + {V2fjg - cos{3(t>)fl)tlEl - 2 sm{3<P)UptgEpEg], (19) 
P'/"' = -^[tlE's-if^-mEl]. 

The SH fields generated by polarization are decomposed into s- and j»-polarized compo- 
nents, 

E^^'' = AgQLeffP^^'' 



and 



E^^^ = Ap^L^fflFsP'^^'' - F^P^^^] (20) 

where VL = 2uj/c is the magnitude of the wave vetor of the SH hght. Ag and Ap (given in 
Ref. [^) are independent of the angle of rotation, but are dependent on the incident angle 
and the optical frequency through a change in the index of refraction. The SH intensity is 
propotional to the absolute square of E{2u!), which for either s- or p-polarized pump beam, 
from (100), (110) and (111) face crystals are found from Table I (only the bulk susceptibility 
is included). 

TABLE I. The second-harmonic field for bulk contibution 
(100) face 

Ep^s = -2d n Leff fc fs tl cos(20)^E2 

Ep^p = dnL,fff,tli2fsF, - /,F,)cos(20)ApE2 

Es,s = 

Es,p = dilFsLeff tlcos{2(t))ApEl 

(110) face 

Ep,, = \dnL,ff [(3cos(20) - l)/2 - 2/2] cos(0)A,E2 

Ep^p = dQL,ff[3 cos2 /2 - fjF, + 2/, sin(0) A^E^ 

Eg^s = SdQ cos sin^ ^sEs 

Es,p = \dQL,ffF,[sm<f)-3sm{3(j))]ApE^^ 

(111) face 

Ep,, = ^dnL,fff^tlsm{3<f))AsE'p 

Ep^p =^,dQ L,ff[V2 f^F, cos(30) + - 2/, - 2f,F,]tlApE^^ 
Es,s = -^dni.ff t-ism{3cj))AsE^^ 

Es,p = ^dQL,ff[V2F,cosi3(P) - F,]tlApEl 

The SH intensity is proportional to the absolute squre of E{2uj). Figure 2 shows 
the rot at ion- angle dependence of the SH intensity Pj„ — Sout for singular ZnSe(lOO) and 
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ZnSe(lll) with only the bulk suscetptibility. For ZnSe(lOO), the SH intensity distribution 
shows a fourfold symmetry, but for ZnSe(lll) a full sixfold symmetry. 



C. Surface or Inrerface Contribution 



The surface or interface nonlinear susceptibility xlurf ^ third rank tensor with 27 
elements, which reflect the symetry of the interface. For SHG the last two indices may 
be permuted at will, thus, the nonlinear susceptibility tensor may have a maximum of 18 
independent nonvanishing elements( e.g. Ci), and higher symmetries lead to a reduction 
in the number of independent and nonvanishing tensor elements. Table II summarizes the 
results for the form of the surface nonlinear susceptibility tensor xlurf various symmetry 
classes. 

(2) 

TABLE II. Independent nonvanishing elements of Xsurf continuous 
point group for a surface in the x — y plane. Where mirror planes exist, 
one of them is assumed to lie perpendicular to the y-axis. [|2ll , p2| 

Point Group Nonvanishing independent tensor elements 

Ci — 1 xxx,xyy,xzz,xyz=xzy,xzx=xxz,xxy=xyx,yxx,yyy,yzz,yyz=yzy, 



yzx=yxz, yxy=yyx, zxx, zyy, zzz, zyz=zzy, zzx=zxz, zxy=zyx 



Cs — m 



xxx, xyy, xzz, xzx=xxz, yyz=yzy, yxy=yyx, zxx, zyy, zzz, zxz=zzx 



C2-2 



xyz=xzy, xxz=xzx, yxz=yzx, yzy=yyz, zxx, zyy, zzz, zxy=zyx 



C2V — Tnm2 



xxz=xzx, yyz=yzy, zxx, zyy, zzz 



C3-3 



xxx=-xyy=-yxy=-yyx,yyy=-yxx=-xxy=-xyx. 



xxz=xzx=yyz=yzy, zxx=zyy, xyz=xzy=-yxz=-yzx, zzz 



- 3m 



xxx=-xyy=-yxy=-yyx, xzx=xxz=yzy=yyz, zxx=zyy, zzz 



xxz=xzx=yyz=yzy, zxx=zyy, xyz=xzy=-yxz=-yzx, zzz 



C^y — Amm 



xxz=zxx=yzy=yyz, zxx=zyy, zzz 



Cq^ — 6mm 



xxz=zxx=yzy=yyz, zxx=zyy, zzz 
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We use the symble x^^J to denote the surface(dipolar) second-order susceptibihty, and 
Px^yfz for surface polarization, e.g., 

TABLE III. Rotation-angle dependences of the second-order susceptibility tensor 
components for different symmetry. 



Surface 


Tensor elements 


P P 

^ m out 


P 

m '-'out 


P 

'-'m out 


'-'m '-'out 


symmetry (100) face xyz 


cos(20) 


cos(20) 


cos(20) 







(111) face xyz " 


cos(30) 


sin(30) 


cos(30) 


sin(30) 




xxz—yyz,zzz 


isotropic 


— 





— 




zxx—zyy 


isotropic 


— 


isotropic 


— 




XXX 


cos(30) 


sin(30) 


cos(30) 


sin(30) 




XXX 


cos(30) 


sin (30) 


sin (30) 


sin (30) 




yyy 


sin(30),sin(/) 


COS0 


cos(30),sin0 


cos(30),cos0 




zxx+zyy 


isotropic 


— 


isotropic 


— 




xxz+yyz 


isotropic 


— 


— 


— 




xxz-yyz 


sin(20) 


sin(20) 


— 


— 




zxx-zyy 


sin (20) 


— 


sin(20) 


— 




zzz 


isotropic 


— 


— 


— 




xxx,xyy,yyx 


cos(30),cos0 


sin(30),sin0 


cos(30),cos0 


sin(30),sin0 




xxz,yyz 


cos(20) 


sin(20) 








xzz 


COS0 


sin0 








xxx,xyy,yyx 


cos(30),cos0 


sin(30),sin0 


cos(30),cos0 


sin(30),sin0 




yyy,xyy,yyx 


sin(30),sin0 


cos(30),cos0 


sin(30),sin0 


cos(30),cos0 




xyz,yxz 


sin (20) 


cos(20) 








zxy 


sin (20) 




sin(20) 






xzz 


sin0 


sin0 







" Only the bulk susceptibility is included. 
The isotropic terms are not shown. 
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= xTj/EyE,, etc., (21) 

in the (x, y, z) system. The transformation rule Eq. (|T3|) is used to transfer the tensor 
Xxy^'^ from crystal coordinates to beam coordinates Xs^I'^- Following Sipe [|l8l, the SH field 
induced by the sheet of polarization can be written as 

(22) 

E;^'f = ApQ[F,e{2uj)Pr^ - F^P^^ 

where e{2uj) is the dielectric constant at frequency 2uj. 

We assume for simplicity that the surface has a simple unreconstructed structure, and 
thus for a particular face it has the same symmetry as the bulk. Therefore we use C^y 



symmetry for the (111) surface and C^y symmetry for the (100) face [|T8[. In fact, there exist 
many local microscopic structures, even in a nominally single surfaces, and in such cases, the 
microscopic symmetry of the structures, such as monoatomic steps, etc. would allow some 
specific tensor element to exist. That is, we may observe a lower symmetry than expected 
if the newly introduced tensor is of different symmetry. Table III summarizes the results in 
cases where surface susceptibility is presented. 

From Table III, we see that there are always plural tensor elements which cause the 
same rotational symmetry. Therefor, the determination of the relevant tensor elements 
must depend on other sources of information, which may be obtained experimentally. Note 
that because SHG is described by a third rank tensor, all surface symmetries higher than 
yield an isotropic response. 



III. EXPERIMENTAL METHODS 

The experimental setup for SHG measurement is shown in Fig. 3. The light source 
of the fundamental frequency is a picosecond mode-locked Nd:YAG (yttrium aluminum 
garnet) laser (Light Conversion Ltd. EKSPLA 2143A) with light pulses of wavelength 
1.06/im, duration time 20 ps, and repetion rate 10 Hz. The linear polarization of the 
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fundamental beam is rotated to the desired angle by the combination of the half wave plate 
(A/2) and the input Glan- Taylor polarizer prism PI. The polarized beam is focused into a 
spot of 3 mm diameter on the sample surface at an incident angle of 45°. In order to avoid 
damaging the sample, the pulse energy was set to about 0.5-0.3mJ per pulse throughout the 
experiment. The visible-cut filter Fl situated immediately before the sample surface removes 
the visible light component, especailly SH signals which are generated in the preceding 
polarizer and lens. The reflected fundamental beam and the SH signal from the surfaceare 
are collinear beyond the sample surface, and the former is eliminated from the beam by 
an IR-cut filter F2 and an interference filter F3 (Melles Griot 03F1B008, A = 550 nm, 
FWHM=70 nm). The linear polarization of the SH signal to be measured by the detection 
system was selected by rotation of the analyzer P2. The signal was finally detected by 
a photomultiplier (Hamamatsu R1464). About 4% of the reflected laser beam is directed 
through the reference arm containing a crystalline quartz plate. The SH intensity from 
quartz plate can be used for providing a reference to remove laser intensity fluctuations. This 
signal is detected by another photomultiplier (Hamamatsu R928) through an interference 
filter F3. Both electronics signals from the two PMT outputs were fed into a BOXCAR 
(SRS 250) averaged over 50 pulses and then recorded using a computer. 

The ZnSe(lOO) and ZnSe(lll) crystals grown by the Markov method were commericailly 
purchased from Ma Teck Material Technologic & Kristalle GmbH, Germany ^3). Their 
surfaces are chemical and mechanically polished. The samples were mounted on a rotation 
stage, with the surface normal set parallel to the rotating axis of the stage and could be 
rotated freely. 

IV. RESULTS AND DISCUSSION 

In Figure 4(a)-(d), the intensity from ZnSe(lOO) is shown as a function of the azimuthal 
orientation for the four different polarization combinations: Pm — Pout, Pin — Sout, Sm — 
Pout, Sin— Sout- AH experiments were made under the same conditions which makes it possible 
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to compare the intensities of the different combinations. Fig. 4(a) [Fig. 4(b)] displays the 
thep- [s-]polarized SH intensity Ipp [Isp] with p-[s-] polarizedincident hght, whereas Fig. 4(c) 
[Fig. 4(d)] shows the s-polarized SH intensity Ips [Iss] with p-[s-] polarized incident light. 
All plots clearly exhibit twofold symmetry. It seems to be consistent with a Cg symmetry. A 
more detailed analysis of the experimental reuslts (see below) shows the interference of the 
surface Cs symmetry with the underlying C^v symmetry of the ZnSe(lOO). The data for the 
four cases are fit by the trigonometric function square plus a constant, this small constant 
can be due to light at other frequencies leaking through the filters. We might be able to 
obtain relative values for the surface and bulk contributions if the function of the incident 
angle multiplying the different surface and bulk contributions would change significantly 
with 6 (see Eq. (^) and TABLE I ). However, with the experimental conditions used here, 
changing from normal to grazing incidence decreases fc and by 5%, whereas fg and Fg 
remain small. Hence any attempt to seperate bulk and surface contribution by the use of 
this method seams to be inappropriate. It is for this reason that in all experiments a fixed 
angle of incidence of 45° was used. 

For simplicity we concentrate on the s-polarized SHG response for an s-polarized pump 
beam, becuase this signal only contains anisotropic terms and will be most sensitive for 
the surface symmetry. This conclusion arises because the bulk contribution is forbidden, as 



seen in Table HI. According to |jT8|, for the (100) face, we see that the surface in general 
has C4t, symmetry. However, such a mechanism should lead no 0-angle dependence. Even 
though we replace the symmetry with lower symmetry C2v, we cannot obtain any 0-angle 
dependent signal. Thus we must allow the effective symmetry of the surface to lower from 
C2v to Cs, corresponding to the observed symmetry of the rotation angle dependence. From 
Table III, for Cg symmetry only the following components of surface or interface nonlinear 
susceptibility do not vanish 

Xxxx, \xyyi tiud Xyxy Xyyx- (^^) 

And a similar experssion for the Sm — Sout case is derived as 
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Es,p{2u;)/El{uj)Ap cc Xxxx sin^ - {xxyy + 2 Xyyx ) sin (j) cos^ </> . (24) 

The solid line in Fig. 4(d) is a least-squares fit of SHG intensity to the data, from which we 
obtain the ratio 

= 0.276 , (25) 

~X.xyy + 2 Xyyx 

the agreement with experimental reuslt justifies our choice of Cg symmetry. There is a large 
scatter in the data for Sin — Sout polarization combination, because its absolute SH intensity 
is small, thus we cannot obtain a complete fit to every data point in this configuration. 



In Ref. dJ, value of Xxxx / iXxyy + 2 Xyyx) — 1-10 was obtained from the NaCl(100)-air 
interface. This difference is not surprising given the fact that the absorption depth in the 
oxide layers of the NaCl(lOO) is more that of ZnSe(lOO) in our experiment. 

This Cg symmetry may be caused by (1) the interface between oxide absorbate and 
sample []7|,|2^, (2) a nonideal cut of the crystal, the domain formation with different symmetry 



axis [|TT|. When one only consider the bulk susceptibility, (f) dependence for Sin — Sout 
is very sensitive to the miscut angle of the crystal axis, even its absolute SH intensity 
is very small. If we consider that the surface tensor elements Xxxx, Xxyy, and Xyxy have 
relevance to the miscut angle of ZnSe(lOO) crystal, which is less than 0.5° PBI], the response 



of Sin — Sout is determined to be more than an order of magnitude smaller than other 
corresponding cases [l^J^. C. Yamada and K. Kimura |^ found it was very difficult to 
measure Sin — Sout rotational- angle dependence because the signal was poor under UHV 
condition. However, we observed that the output was about half of the other cases. So we 
concluded that in our experiment the interface plays an impoartant role for producing the 
anisotropic in Sin — Sout configuration. This is supported in part by similar measurements 
made for GaAs(lOO) surfaces contaminated with carbon and oxygen which revealed a high 
degree of anisotropy in the SHG signal The anisotropiy cannot be ascribed to any 

particular surface symmetry. One possible explanation may be suggested: the anisotropy 
may arise due to variations in oxide layer thickness which could be regarded as generating 
a stepped buried ZnSe( 100) /oxide layer interface. 
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For three configurations {Pin — Pout, Sin — Pout, Pin — Sout), the SHG intensity distribution 
shows a 6*4^, symmetry if we only consider the bulk susceptibility. Further, we notice that 
the oxide layer interface would lower the effective symmetry of the surface from C^^ to 
C2V In Fig 4.(a)-(c), we observe twofold symmetries not fourfold for three configurations, 
which is the result of interference of bulk SHG with surface SHG. Threfore, the observed 
rotation-angle dependence can be described as 

Is.p(x\A + B cos(20) + C sin(20) f , (26) 

where A represents the isotropic surface contribution, B is the bulk dipole term and C the 
surface anisotropic term. We least-square fitted the observed rotation-angle dependence for 



three configurations by Eq. (p6D with A, B and C used as adjustable parameters. These 
fitting parameters are summarized in Table IV for the SHG data. As shown in Fig. 4(a)-(c), 
the fit was satisfactory. 

TABLE IV. Values obtained for fitting parameters in Eq. (p6D 



Configuration ABC 

Pin -Pout -0.034 0.203 1.201 

Pin-Sout 0.046 1.200 -0.126 

Sin- Pout 0.062 0.039 1.150 

In the Pin — Pout combination (Fig 5.), the peak at = 135° is higher than that at 
(f) = 45°, and the peak at = 315° is higher than that at = 225°. These diffrences arise 
from the interference between the isotropic surface SHG and the bulk SHG. The signs of 
the nonlinear susceptibility components for the exciting electric fields in the (Oil) and (Oil) 
planes are opposite. This is because the bulk Zn-Se-Zn-Se- chain along the [Oil] has Zn 
atoms higher than Se atoms and Zn-Se-Zn-Se- chain alone [Oil] has Se atoms higher than 
Zn atoms. This phase difference of 180° between the electronic wave functions of the two 
bulk chains leads to the diffrence in the interference between the bulk and surface SHG and 
to the change of the SHG peak intensities. 
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V. SUMMARY 



We have developed a phenonmenological model and performed experiments to determine 
the symmetry of the noncentrosymmetric semiconductor crystal ZnSe(lOO) covered with an 
oxide layer. Under the four combinations of fundamental and harmonic linear polarization 
states considered, the harmnic intensity can be experessed as a function of the light polar- 
ization and as a function of the azimuthal angle. Using the caclulated sensitivity to rotation 
angle of the SHG signal from the bulk and surfaces, we have been able to deduce the sym- 
metry of the noncentrosymmetric crystal surfaces. We found that the measurement using 
Sin — Sout is particularly useful in determining the symmetry of the oxdized layer interface, 
which would lower the effective symmetry of the surface from C^y to C2v The separation be- 
tween bulk and interface or surface SHG demonstrated here is promising for the application 
of this technique to the study of surface and interfacial properties. This way can be used to 
detect the quality of the substrate surface for growth of II- VI epilayers on ZnSe layers. 

In addition, we have shown that the [Oil] and [Oil] directions can be distinguished 
through the analysis of p-incident and p-output confugration. 
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FIGURES 

FIG. 1. Principal geometry for SHG in reflection. In the experiment the surface is rotated its 
normal Z (azimuthal angle (p )> and the polarization of both fundamental and SH can be varied to 
any direction between s and p. 

FIG. 2. The rotation-angle dependence of the SH intensity Pj„ — Sput for singular ZnSe(lOO) 
and ZnSe(lll) with only the bulk suscetptibility. 

FIG. 3. Schematic diagram of the experimental configuration used to measure the SHG from 
ZnSe single crystal surfaces. 

FIG. 4. Rotational anisotropy of SHG from ZnSe(lOO) surface in free air for a set of input and 
output polarization combinations. The circles are the experimental points, and the solid lines are 
least-squares fit to theoretical calculation. 

FIG. 5. Polar plot of the SHG intensity from ZnSe(lOO). The pump and SHG beam were both 
p polarized. 
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